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Resume 

Dans ce travail on a propose et etudie une formulation mixte a trois champs pour resoudre le probleme de Stokes avec 
des conditions aux limites non-lineaires, du type Tresca. Deux multiplicateurs de Lagrange ont ete utilises afin d'imposer 

div{u) = ct dc rcgulariscr la fonctionncUc cncrgic. Lcs elements finis PI bulle/Pl-Pl ont pcrmis de discretiser le 
problcmc resultant. Dcs estimations d'crrcurs ont etc dcrivccs ct plusicurs tests numcriqucs sont realises. 

Mots cles; Probleme de Stokes, Frottement de Tresca, inequation variationnelle, elements finis mixtes, estimation 
d'erreur. 

Abstract 

In this work we propose and study a three field mixed formulation for solving the Stokes problem with Trosca-type 
non-linear boundary conditions. Two Lagrange multipliers are used to enforce div{u) = constraint and to regularize 
the energy functional. The resulting problem is discretised using PI bubble/Pl-Pl finite elements. Error estimates are 
derived and several numerical studies are achieved. 

Key words/ Stokes problem, Tresca friction, variational inequality, mixed finite element, error estimates. 



Notations 

We need to set some notations and recall some functional tools necessary for our analysis. Let Q C R'', d > 2, an open 

set with boundary dQ wich is the union of two nonoverlapping portions To and F (may be empty). 

The euclidian norm of a point x £ R'* is denoted by |x| in what follows. The Lebesgue space L^(f2) is endowed with the 

norm: 

WpeL^Q) ||p||o= (^^|p(x)|2dx 
while Lq{Q,) is the closed subspace of L^(f2) defined by: 

Lo(a) = |p e L^{Q) such that J p(x) dx = o| . 

We make constant use of the standard Soblev space H"^{fl), m > 1, provided with the norm: 



/ \ 1/2 

E I 

I 0<|a|<m 



mil 



where a is a multi-index. Fractional Sobolev spaces H''(Q), v e R_|_\N are defined by 

= {(^ e such that ||<p||,,n < +00} , 

with 

\ \a\=m I 

with m being the integr part of v and Q its decimal parts. 

The closure in _H"'^(n) of -D(n), the space of infinitely differentiable functions with support in f2, is denoted by Hq{Q). 
On any portion F C dO, we introduce the space H 2 (F) as follows 



iJ-5(F) = |(p e L^(F) such that Hi^Hi^r < +oo| , 



where 

^2 



The space H 2 (r) is the dual space of (F), (•, •) stands for the duality pairing and 

M-hr= sup 

1 1 
The special space iJQQ(F) is defined as the set of the restriction to F of the functions oi [dO.) that vanish on 9f2\F 

1 

and its dual space is denoted by (i?Qo(F))'. 

The cartesian product of k previous spaces and their elements are denoted by bold caracter. The respective norms 
are introduced as follows: 

N 1/2 



i=0 / 

k \ 

Ell/^ill-ir J M= (mi,--- € H-5(F), 

i=0 / 



Let X C H^{Q.) be a subspace of functions vanishing on a non-empty portion Fq open in dO, 

X =^ve H^{Q) such that v^^^ = o| . 
We introduce the enrgetic norm 1 1 1 • 1 1 1 1 in X corresponding to the scalar product 



V £ij(u)£ij(v)d(x), 
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where Sij is the ij-th component of the hnearized strain rate tensor e(u) = g V*u). From the Korn inequahty it 

follows that II ■ 111 and ||| • |||i are equivalent in X. 

We denote by n the outward unit normal to d^l and u„, respectively ut, the normal , respectively the tangential, 
component of u. 

The stress vector a is equal to a.n where a is the Cauchy stress tensor defined by: 

a — 2i/e{u) — pS, 

where p is the hydrostatic pressure, 5 is the identity tensor and u is the kinematic fluid viscosity. 



1 Introduction 

No-slip hypothesis at fiuid-wall interface leads to good agreement with experimental observations for newtonian fluids 
which is no longer true for non- newtonian fluid T. For example, in the flow of certain high molecular weight linear 
polymers through circular dies, the exit flow rate has been found to be a discontinous function of pressure drop over a 
certain range of shear rates [21 [3]. This obervation is consistant with the hypothesis that the velocity at the wall is not 
zero. Several studies have been made and showed not only that slip takes place when a threshold is reached [1] but also 
it's the origin of many defects and instabilities in the polymer injection process [5l[6]. 

The first attempt to integrate this boundary condition in a numerical simulation of a fiow is due to Doltsini et al. 
[7] and Fortin [8j. Since that, many papers were published simulating various fiows with such boundary conditions (see 
[9] and refrences therin). Recently, based on the penality method, error estimates for the Stokes problem with Tresca 
boundary conditions with strong regularity assumption on the velocity field are obtained |10j . 

The aim of this work is to contribute to the numerical analysis of Stokes problem with Tresca boundary conditions. 
Our first purose is to carry out the convergence analysis and a priori estimates for the mixed finite element formulation 
of the above cited problem. The second one is to derive an algorithm well adapted to this formulation and easy to 
implement in order to validate our theoritical estimates. 

The paper is organized as follows. First, we introduice the equations modelling the Stokes problem. Then we 
establish the continous mixed variational formulation is section 3. The following section is devoted to a priori error 
estimates , we show an optimal order of h?'^'^ with H^(r2) assupmtion regularity on the velocity. In section 5 we propose 
an algorithm based on augmented lagrangian method to solve the 2D problem and make some numerical tests. 

2 Setting Stokes problem with nonlinear boundary conditions 

We consider the following Stokes problem with nonlinear boundary condition of Tresca friction type: 

—div(ve{\i)) -I- Vp = f in SI 

div{\i) = in f7 

u = on Po , , 

u„ = on r ^ ' 

\<Jt\ < g ^ Ut = on r 

\at \ = g ^ 3fc > a constant such that Ut = — fccrt on F 

with n C K'* (d = 2 or 3) an open set with regular boundary dO,, which is the union of two nonoverlapping portions Fq 
and F. Fo is subjected to no-slip boundary condition while F is where le fluid may slip. We need this result to derive 
the variational problem. 



Proposition 2.1 \1_ 



^ , , . ,. , , (Tt.ut -I- g|ut| = 071 F (2) 

at \ — g ^ Mt ~ —Kat k IS a non-negative constant on i 

One can derive the variational formulation of ([T]): 

{Find u G Vdiv{^) such that :Vv G Vdi„(n) 
(3) 
a(u, V - u) + j(v) - j(u) > L(v - u), 

with 

V{Q.) = {v G H'(f^), V|ro = 0,v.n|r = 0} V^^n) = {v G V(n) , div{v) = in fi} 

a(u,v)= / !/e(u) : e(v) dr2 L(v) = / fvdfi g a non-negative function in I/^(F) 
Jn Jn 

and = J g\^t\dr Vu,v G Vd,„(n). 

Problem ^ is an elliptic variationnal inequality of the second kind which has a unique solution T3j. Moreover, since the 
bilinear form a(-, •) is symmetic ([3| is equivalent to the following constrained non-differentiable minimization problem: 
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Find u £ 'Vdiv (f^) such that : 

(4) 

J{u)<J(v) VvGVd„(n), 



where J{v) = ^ a(v, v) + j(v) - L{v). 



3 Mixed Formulation 



In order to solve Q a Lagrange multipher q is needed to enforce the condition div{u) — in which can be identified 
with the pressure. In the other hand Fujita proved in [19 that (|3| is equivalent to 

f 3(7t G (H(|)(r))' such that \at\ < 5 on F 

(7t(v-u)t+j(v)-j(u) > VvGV(f7). 



at is seen as a Lagrange multiplier and can be identified with the shear stress on F. The minimization problem Q is 
equivalent to the following saddle-point formulation : 

(Find (u,p, X) eH such that: 
(6) 
■C(u,g,^) < £{u,p,X) < £(v,p, A) V(v,g,At) G H. 

C{v,q,^i) = la{y,y)- I q div{v) + [ fj.Vt - L{v) V(v, g, ^) G H = V X Lg(f7) X Q 
^ Jn Jr 

where 

Q = € (L'(F))'*-M/il < 3} 

^ie{L\^))^-\ ^/ivt - ^3|vt| < ovv g v [ 

According to |20| problem ([gJi has a uniqe solution charcterized by 

Find (u, (p, A)) G V X A such that: 

a(u,v) + 6((p,A),v) = L(v) Vv G V (7) 

[ 6((g-p,M- A),u) < V(g,M)eA, 

where 

b{{p, A), v) = -(p, dj?;v) + (A, vt) (8) 
and A = Lg(n) x Q is a closed convex of X = Lo(!^) x (^^(F))''-^ 

Lemma 3.1 There exists a constant a > such that : 'i{q, fi) G A4 

snp'-%^>a(M + MA. (9) 

v6V l|V||l,Q V 2/ 

Proof : To prove this result we are inspired by ^0]. We have to prove that for all {q,n) G A4 there exists u G V 
such that: 

div u — q in Q 

u = onFo 

u„ = on F ^ ' 

ut = h^^{fi) on F 

wich satisfies 

llulli < C (||g||(, + ||At||_i) , (11) 

-1 - - 

where h (■) is the inverse Riesz operator (HgQ(F))' — >■ HQg(F). 

The proof is diveded into five steps: 

Step 1 
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We suppose that J7 is a convex with a regular boudary dfl. Let q G Lo{Q,) and <j>i be the solution of: 



q on SI 
in do. 



According to 1211 problem (12 1 admits a unique solution 4>i verifing 4>i G H^{Q,) and 

ll0i||2 < C||g|lo 

where C is a constant independent of and q. 
Step 2 

Since <l>i e H^{n), ^ e H^idQ.), we now consider the following Neumann problem: 



A02 

d(t)2 
dn 



dn 



on f2 
in do. 



According to pT], this problem admits a regular solution 02 G verifing: 

IWbSCIlf II.. 

step 3 

Let ^ be the unique solution to the following bilaplacian problem : 



dn 



= in S7 
= in an 

= X in 



where 



X = 



«^ "'^'^^0 

or OT 

d(j)i 902 , ,-1. . ^ 
or OT 



we can easily show that x ^ (d^). From [22] it holds 
Step 4 

Setting u = V0i + V02 + cur/ ip with 



lV^ll2<cilxlU. 



dip 

dX2 

dtp 
dxi 



so that ut = h ^{fJ.) on F. 
Furthermore we obtain: 



l|ul|i < ||V0i||i + ||V02||i + ||curZi/)jji 

< ||0l||2 + II02II2 + llV-lb 

< cfiMio + II^IU + llxlU 

V an 2 "2 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



where C > is a gene ri c co n stan t. 

Using in equ alities ( 13 1, (jlSl, (|l7|, the continuity of normal trace application from H^(f2) onto Hs (9f2) and the continuity 



of ft. ^, (|l9|) becomes 
Step 5 



u||i < C ( IIqIIo + IImII- J 
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sup ^(('?'^)''^) > ''((g^/^)'") 

vGV ||vl|i,n - ||u||i 



(20) 

1 IMl' + llMiP_i 
- C \\q\\ + M\_i' 

> ^{M\ + \ML. 

Then take a — -\- to finish the proof. □ 
2C 

Theorem 3.2 \20l Suppose that a{-, ■) is continuous, V-elliptic bilinear form on V(r2) and ^ holds. Then there exists 
a unique (u, (p, A)) solution of mixed problem Moreover, (u, (p, A)) is also the unique solution of the saddle-point 
problem 

4 Error estimates 

The present section is devoted to finite element approximation of the saddle-point problem ([6|. The key point lies in 
finite element discretization of the closed convex Q of the Lagrange multipliers which leads to a well-posed discrete 
problem and gives a good convergence rate for the approximate solution. 

We use classical PI bubble-Pl finite element to disretize (u,p) and PI finite element on T for the Lagrange multiplier A. 
This choice is motivated by Brezzi's and Sassi's results, see [291 130| . 

O is supposed to be polygonal. Let 7h be a regular partition of SI with triangles in the sense of [24]. We denote by Pn(K) 
the space of polynomials of degree less and equal to n G N defined on k £ Th- We denote by Bk. the space of bubble 
functions defined on k which is a sub-space of Hq{k). Then we can define the following discrete spaces : 

B = Vh = {vh G C°(n);vh\^ G Pi Vk G Th, Vh\ro = 0, and v,,.n|r = 0}, 

Vh = [Vh + B]^ Wh = {Vfcir, Vh G Vh}, 
U = {qh G C°(n); gh|« G Pi Vk G Th, I qn = 0}, 

Qh = l^ft G Wh, I t^h^h - I g\M < OVV'h G Whj , 
Alh = Lh X Wh, Ah = Lh X Qh. 

Remark 

Qh is an external approximation of Q, so the discretization is non-conforming and would weaken its convergence order. 
Disicretizing ([7| we obtain 

Find (uh, (ph. Ah)) G Vh x Ahsuch that: 

a(uh,Vh) -hb((ph,Ah),Vh) = L(vh) Vvh G Vh, (21) 

b{{qh — Ph, IJ.h — \h),Uh) < '^{qh, fJ-h) € Ah, 

where Ah is a closed convex of A^h- 

A sufficient condition for the existence and uniqness of the solution to problem ( |21[ | is the inf-sup condition [29j . 

Lemma 4.1 These two propositions are equivalent: 

* There exists a constant /3 > independent of h such that : 

K[lh,tJ-h),^h) 



sup 



,VhJ >/j(^||g^||„ + ||^th|Li) 'i{qh,t^h)eMh. (22) 

Vh 1 V 2/ 
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* There exists two constants /3i > and /32 > independent of h such that 

{qh,div\-h) 



< 



sup ^ rr—^ > Pi\\qh\\o yqh&^h 

v^ev,, ||v;i||i 

and 

\l-^h,'^th) ^ o II II w ^ TA, 



Proof : 

To prove this result it suffices to show that: 



sup {(/ihjVth) + (g;j,di?;(vh))} = sup {(^h,Vth)}+ sup {{qh,div{^h))} 



Let us prove the non trivial direction 



sup Vth)} + sup {{qh,div{vh))} < sup {(/x;,, Vth) + (g^, di?;(vh))} 



which is equivalent to 



(phjUth) + (gh,diw(wh)) < sup {(Ath, Vt^) + (5^, dii;(vh))} V(u,i, w^) G V,, 



(23) 



(24) 



(25) 



We suppose that (251 is not valid, ie 3 {6fi,ujh) G V;^ such that 



sup {{fJ.h,Vth) + {qh,div{vh))} < {fih,Oh) + {qh,div{LOh)) 



(26) 



(261 => {fih,vth) + {qh,div{vth)) < {^h,Oh) + {qh,div{u)h)) Vu/i G Vh 



{^Ji.h, vth — Oh) + {qh, div{vh — i^h j) < 0. 



Since we know that: 



e (Llir))"-' ^\\^^^\< sup i^^^^ ^ 



V^h G {Lliryf-^ 3lp,^ G Vh such that 



(27) 



and suppose that (fih, Vth — Sfe) < Vuh G V;i, wich remains valid for Vh = 'i^h + Oh so that 

{fJ'h,^h) < 0, 

which contradicts ( |27[ ). 

The same reasoning can be applied to {qh, div{vh — i^h)) < since (■, div(-)) verifies similar inf-sup condition. This ends 
the proof. 

□ 



Proposition 4.2 There exists a constant /3 > independent of h such that 



b{{qh,fJ-h),Vh) 



sup > 

^hGVh l|Vfe||l 



/3(l|gh||o + ||Mh||-i) V((7h,/ih) G A4h 



Proof : Recall that 

y {(qh,fih) ,Vh) G Mh X Vh b{{qh,^h) ,Vh) = {-qh,div{-Vh)) + {fih,^th) 

According to lemma [4Tj it suffices to show that : 

=1 fl ^ n 4- 1 w ^ TT (Qh.diWh) 
a /5i > tel que V gji G L/i sup — t-\ n > Pi||gh||o 



and 



3/32>0 telqueV/ihGWh sup > /SallMfell 

vftSV^ ||v;i||i 



which are both ( 29 1 and ( 30 1 established in and [35j respctively. 

Now we will derive error estimates for primal variable, being inspired by [27| . 



(28) 



(29) 

(30) 
□ 
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Lemma 4.3 I26f Let {u,p, A) and {uh,ph, Ah) be solutions to respectively. Then for any (vh, qh, fih) G Vh x Ah 

it holds: 



a(u - Uh, u - Uh) < a(u - Uh, u - Vh) + b{{p, A) - {qh,p.h), Uh - u) + fe((p, A) - {ph, Ah), u - Vh) 
+6((p,A) - ((7h,Mh),u) + 6((ph,Ah) - (p,A),u). 
Proof : Let Vh be an element of Vh. It follows that: 

a(u - Uh, u - Uh) = a(u - Uh, u - Vh) + a(u - Uh, Vh - Uh). 
Using the first equations of ([7| and of (21 1, this gives : 

a(u - Uh, Vh - Uh) = a(u, Vh - Uh) - a(uh, Vh - Uh), 

= L(vh - Uh) - b{{p, A), Vh - Uh) - -L(vh - Uh) + b((ph, Ah), Vh - Uh), 
= A),Uh - Vh) + &((ph. Ah), Vh - Uh). 



(31) 



Then we deduce 



a(u - Uh, u - Uh) = a(u - Uh, u - Vh) + b{{p, A), Uh - Vh) + 6((ph, Ah), Vh - Uh). 



Finally we have 

a(u - Uh, u - Uh) = a(u - Uh, u - Vh) + b((p, A) - {qh,p.h), Uh - u) + fe((p. A) - (ph, Ah), u - Vh) 
+b((p,A) - (gh,Mfe),u) + - (p,A),u) 

+b{{qh,fih) - (Ph, Ah), Uh). 



But according to (21 1, b{(qh, iJ.h) — {ph, Ah), Uh) < for all (qh, ^h) G Ah. This ends the proof of the lemma. □ 
We now derive an upper bound of the terms involved in ( |31[ ). 
Lemma 4.4 Let (u,p. A) and (uh, Ph, Ah) be solutions to \21^ respectively. Suppose that u £ H'^(f2) andp G 77^(0). 



T/ien 



||u- Uhll? < C(u,p,.g) (^/i||A- AhlLi j 
where C{u,p,g) is a positive constant depending only on ||u||2, ||p||i and \\g\\L'^(r)- 



(32) 



(33) 



Proof : Using Lemma 4.3 we will show that there exists (vh, {qh,fJ-h)) G Vh x Ah satisfying : 
a(u — Uh,u — Vh) < C(u)ft||u — Uh||i, 

b{{p,X) - {qh,iJ,h),Uh - u) < C(u,p)ft||u - Uh|ji, 

b((Ph,Ah)-(p,A),u-Vh) < Ch{\\u-Uh\\i + C{p)h + C{u)h} 

b((p,A)-(gh,Mh),u) < C{ufh^ 

b((Ph,Ah)-(p,A),u) < C(u) (ftllA- AhlUi +C(u)/ii +C(g)/ii) 



Before proving these estimates, we first have to recall some useful results. Let Ih, 3h and ih be the Lagrange 
interpolation operators on Vh, Lh and Wh respectively. From [24], there exists a positive constant C such that 
Vv G H2(n), Vp G Lh and Vt/j G i/3(r): 



|v-IhV||i <Cft||v||2, |b-3hP||o < ChMi, W-ilMoT < Ch2\ 

Let TTh be the projection operator from (I/^(r))'^~^ on Wh defined by: 

TThV- G Wh, y (TThV- - -0) = V^hGWh. 

It holds Vr G [0, 1] and G [0, r + |] one has: 

VV- G /i-5|j^_7rhV'||_i r + ftni^-'rhV'IU,r < C/i"+i 



(34) 



(35) 



(36) 



Let Ilh be the projection operator from L^{Q) on hh defined by: 



Finally, let us note the trace theorem implies that 



{nhq-q)sh = VshGLft. 



|A||i,r < C||u||2 



(i) The first term is evaluted by using the continuity of a(-, •) and the property (34 1 



(ii) Using ( 34 1 and ( 38 1 we have 

HiP,^) - {qh,fj.h),iih - u) 



i{u- Uh,u-~ Vh) < Cllu - Uh||i ||u - Vh||i Vvh G Vh, 



< C inf {llu — Vhlli} llu - Ufelli, 



< C||u — Ihu||i||u- u/i||i, 

< C(u)/i||u - u,j||i. 

= -(p - qh, div{uh - u)) + (A - fih, um - ut) , 

< c|llp-gfello + ||A-/i,,||_i|||u-Uh||i '^{qh,Hh) e J^h, 



< C<^ inf \\p~qh\\o+ inf ||A - i M|u - u^iii, 



< C{u,p)h\\u- Uh\ 



(Hi) Further 



&((p, A) - (ph, Ah),U" Vfe) = -{p-ph,div{u-Vh)) + {X- Xh,ut-Vht) , 

< c||b-p,,||o + ||A-Ahll_i j||u-Vh||i VvfeeVfc, 

< c{|b-g;,||o + ||A-/i„|l_i_r 

+ \\qh -Ph\\o + \\Xh - MhlLi,r} II' 



lu - Vhlli. 



Since {ph — qh, Xh — fJ,h) & Mlh then it follows from the discrete inf-sup condition (281 

bjiPh - qh,Xh - Mfc) >Vh) 



13 (j\ph — g;i||o + ||Afe — /ih||_ i pj 



r I sup 



1 



< sup 



f b{{ph-p, Afe- A),Vh) _^ b{{p- qh,X- /ih) ,Vh) 



a(u- Uh,Vh) 



— 

vhSVfe l|v;i||i 



+ l|p-gfe!lo + ci|A-ph||_i 



Hence : V (qn, fJ-h) G Alh 

||ph — + ||Afe - /ih||_ i_r < C* ^ + Hp — gh,||() + ||A - i^P^ 
Then combining the last inequality of (391 and (41 1 and using property (34 1 we obtain 

b((p, A) - (ph, Afe), u - Vh) < C<||u-Uh||i+ inf \\p - qh\\o 

+ inf l|A-ph||_il inf ||u-Vhl|i, 

< Ch{\\\x-Uh\\i + C{p)h + C{\x)h} 
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(iv) To estimate this term we invoke the definition of the L^-projection operator: 

fe((p, A) - (gh, ^h), u) = -{p-qh,divu)+ {X-fih)utdr V/ih € Qh, 



= y (A - 7r,iA)utdr for = tt^A, 

= y (-^ ~ 7rhA)utdr ^ j (^^ TTh\)7VhUtdr, 

= J {X - nhX){ut - nhUt)dr, 

< \\X - 7rhA||j,2(r)llut - 7rhUt|1^2(p), 

< C{ufh^. 

(v) Now we shall estimate the fifth term of (33 1 using (|2| and the definition of Qh 



(42) 



r 



Ah) - (p, A), u) = ~{ph - p.divu) + I {Xh - X)vLtdV, 
(Ah - A)utdr, 



r 



= j (Ah - A) (ui - ih(ut)) + I (Ah - A)ih(ui) + j [Xut - g\\it\ 

< J (Ah - A) (ui - ih(ut)) + J g {\ih{ut)\ - \ut\) + J X (ut - ih{ut)) 

< y (Ah - A) (ut - ih(ut)) + J g\ihiut) - ut\ + J X {ut - ih{ut)) , 

< \\Xh - A||_ i_p||uf — ih(uf)|| i_p + ||g||o,r||ut — ■jh(uf)||o,r 
+ ||A||(),r||ut — Zh(ut)||o,r, 

< C(u)/t||A- Ah||_i,r + C(u)2/ii +C{g)C{u)hi, 

< C(u,g)|;i||A- Ah||_i^r + /i5|. 



Assembling the estimates (i)-(v) in the Lemma 4.4 and using the V-ellipticity of the bilinear form a(-,-)i finally 
arrive at the following estimate 

||u - Uh||? < C(u,p) (^h\\X — Xh\\_i_ + /i||u — Uhllj + C(u, g)/i2 , 

then using the Young inequality we can write for every constant /? > 

C(u,p)/i||u - Uh||i < C(u,p) (jSh^ + i||u - Uhll? 



Taking /3 such that < 1 then leads to the desired result. □ 

Lemma 4.5 Let (u,p, A) and (uh, Ph, Ah) he solutions to \21^ respectively. Suppose that u G H'^(n) andp £ H^{Q). 
Then 

\\p-ph\\a + IjA- Ah||_i < C(u,p){/i + ||u-Uh||i}, (43) 
where C(u,p) is a positive constant depending only on ||u||2 and ||p||i. 



Proof : Using (40 1 and (41 1 we get the disired result. □ 



Theorem 4.6 Let (u,p, A) and (uh,Ph,Ah) be solutions to \2l]) respectively. Suppose that u £ H^(f2) and 

pe H\n). Then 

llu-Uhl|i + llp-Phllo + HA- Ahll_i < C{u,p,g)hT, 
where C{u,g) is a where C{\i,p,g) is a positive constant depending only on ||u||2; ||p||i and \\g\\i^2(^Y)- 
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Proof : By assembling ( 32 1 and ( 43 1 we can write: 

||u-u,,ll? < C(u,p,ff) |ft||A- Ah||_i 

< C{u,p,g)h{h + \\u- Uh\\i} + C{u,p,g)hi 

< C{u,p,g)h'^ + C{u,p,g)h\\u - Uh\\i + C{u,p,g)hi 



then using Young's inequality we can easily write: 



so that 1 43 1 becomes 



|u - Uh||i < C{vi,p,g)h* 
b- Phllo + P - Ah||_ 1 < C(u,p,g)/ii 



wich leads to the desired result. 



□ 



5 Numerical simulations 

We briefly describe the numerical resolution of the 2D Stokes problem with boundary conditions of Tresca friction type. 
For this aim, the augmented lagrangian method [31] will be used. 
The minimization problem Q is replaced by : 



Find (u, $) G n such that : 

E(u,$) < E(v,¥?)V(v,^) e n, 



where 



n = {(v, tp) G Vdi„(r2) X L^(r) such that = Vj}, 
and E the lagrangian is defined on 11 by: 

V(<^,v) G n E(v,<^) = ^a(v,v) -L(v) + 

Then, the following saddle-point problem is derived 

{Find (u, A) G n X L^{V) such that: 
Cr(vi,^,^l) < £.(u,$. A) < Cr{v,<i>,X) V(v,v9,m) G Vd„(f7) X {L^(r)f, 

where 



(44) 



(45) 



(46) 



and we use bloc relaxation Uzawa algorithm, or ALG2 as mentionned in |31| , to solve (45 I. This leads to the following 
algorithm: 
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1. Initialisation: il> \ \° et r > fixed. 

2 . Repeat until convergence 

Find u'° G V such: Vv G V 
a(u\v) + r(u^vOr = L{^r) + (r^"-^ - X\vt) 
div{u'') = 



(47) 



\\^+ru1\\ - 
rllA*^ + rUf\ 



if 



unless 



l|A^- + ru?l|>5 



(48) 



\k + l 



(49) 




Remarks 



* It's recommanded in [311 132| to choose pk = p = r to ensure the convergnece of the above algorithm; 

* A second issue is how to choose r? Numerical tests show that there is an optimal value Vopt for which convergence 
is the fastest. Unfortunately, this result still unprooved. 



5.1 Numerical Tests 

A no-slip 2D Stokes solver [33] is used and Tresca friction boundary conditions were implemented on. Q is the square 
[0, 0.1]'^, the fluid can slip on F = Tupper U Viower = [0, 0.1] x {0.1} U [0, 0.1] x {0}, the viscosity is taken equal to 0.1 and 
10~^ is choosen as a stopping criterion. 



5.1.1 Test 1: 

If the threshold is never beeing reached then there is no-slip on all parts of the boundary dn wich is the case if the 
solution (u, p) is that of the Stokes problem with homogeneous Dirichlet boundary conditions. Logically, the value of g 
has no effect on the solution of such problem. 
The volume data f is adjusted so that the exact solution will be : 

ui(x,y) — —cos{20nx)sin{20Tvy) + sin{20TT'y) 
U2{x,y) — —sin{20nx)cos{20Tvy)~sin{20TTy) 
p{x,y) — 20n{cos{20-Ky) — cos{20-Kx)) 

As shown in table [l] error decreases as we consider smaller mesh size. 
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l/h 



np 



nt 



|u - U/j||o 



u 



U/illl \\P-Ph\\0 



700 
900 
1100 
1300 
1500 
1600 
1700 
1800 
2000 
3000 



8522 

14038 

20880 

29506 

39103 

44756 

50228 

56385 

69068 

155610 



16762 

27714 

41318 

58490 

77604 

88870 

99774 

112048 

137334 

310018 



1.245e-04 
7.446e-05 
4.969e-05 
3.534e-05 
2.647e-05 
2.306e-05 
2.062e-05 
1.837e-05 
1.514e-05 
6.650e-06 



1.365e-01 
1.056e-01 
8.656e-02 
7.305e-02 
6.304e-02 
5.889e-02 
5.554e-02 
5.259e-02 
4.762e-02 
3.159e-02 



4.253e-03 
2.992e-03 
2.211e-03 
1.700e-03 
1.451e-03 
1.317e-03 
1.262e-03 
1.170e-03 
9.956e-04 
6.181e-04 



Table 1: h: mesh size, np: number of noeuds, nt: number of triangles 
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l|u - u/illo 


\W - u/illi 


\\p-Ph\\o 


nit 





3.0405e-03 


9.1986e-02 


7.4967e-02 


26 


0.015 


3.0273e-03 


9.1623e-02 


7.7007e-02 


131 


10 


3.0251e-03 


9.1596e-02 


7.7089e-02 


135 


40 


3.0251e-03 


9.1596e-02 


7.7089e-02 


135 



Table 2: Effect of g on the approximate solution, r = 10, nj^inumber of iteration to convergence 



5.1.2 Test 2: 

We set g — 0.015 wich is consistent with experimental values, see [3] and [11], and we enforce parabolic profil on both 
Fje/t and T,.ight'- 



Ul =Ur = 



y(i - y) 
- y) 



where u; = u|r,,^.j and = u\r^^^^,. 

We choose this profile to enforce shear stress near the solid wall to reach the threshold without considering a complicated 
domain geometry. We can easily notice that fluid slips on some regions of dQ and adheres the other regions, see figures 
(ITl|2l. 
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Figure 2: Zoom of snon-slip and slip zones 

Since an explicit solution to such a problem is not available, we calculate the discrete solution with sufficiently refined 
mesh, h = jij^i which is taken as the reference solution; next we compute u^, the approximate solution, for different 
mesh sizes h and we compare them to the reference solution. 



h 




u - UhWo 




ao 




U - Uh\\l 




ai 




p-Ph 







Op 


1 


4286e-03 


3 


5378e-04 


1 


213 


8 


8945e-03 





720 


2 


0119e- 


01 





244 


1 


2500e-03 


2 


7779e-04 


1 


225 


7 


6057e-03 





729 


1 


7985e- 


01 





256 


1 


lllle-03 


2 


2973e-04 


1 


231 


6 


8531e-03 





732 


1 


6176e- 


01 





267 




l.e-03 


1 


8160e-04 


1 


247 


5 


9045e-03 





742 


1 


4513e- 


01 





279 


9 


0909e-04 


1 


5163e-04 


1 


255 


5 


3128e-03 





747 


1 


3040e- 


01 





290 


8 


3333e-04 


1 


2746e-04 


1 


264 


5 


0581e-03 





745 


1 


1660e- 


01 





303 


7 


6923e-04 


1 


0930e-04 


1 


272 


4 


6551e-03 





748 


1 


0512e- 


01 





314 


7 


1429e-04 


9 


5356e-05 


1 


278 


4 


3974e-03 





749 


6 


0087e- 


02 





388 


6 


6667e-04 


8 


5856e-05 


1 


280 


4 


6691e-03 





733 


8 


9121e- 


02 





330 


6 


2500e-04 


7 


3662e-05 


1 


289 


3 


7694e-03 





756 


3 


9288e- 


02 





438 


5 


8824e-04 


6 


5507e-05 


1 


295 


4 


1820e-03 





736 


5 


8219e- 


02 





382 


5 


5556e-04 


5 


8187e-05 


1 


301 


3 


6802e-03 





747 


4 


9130e- 


02 





402 



Table 3: Convergence rates with respect to h 

Table [sl provides the variation of ||u — Uh||o, ||u — Uh||i and ||p — pfi||o with respect to the mesh size respectively. 

3 

The first remark one can make is the rate convergence of _ff^-norm of error on u is equal to - which is in agreement with 
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theoretical result. The second one is that in spite of concidering very small mesh size, h = j^^i we cannot conclude 
about rate convergence of u and p error L^-norms. 



6 Conclusion 

A three field mixed formulation of the stokes problem with Tresca boundary condition has been introduced and studied. 
The convergence analysis and a priori error estimates of the discrete corresponding problem have been established. In 

3 

particular, wc show an optimal error estimate of order /i4 for the velocity when it is approximated by classical PI bubble 
finite element. A numerical realisation of a model example have been proposed wich confirms the theoritical result. 
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